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Abstract 

We study the configuration space of the Tomboulis £0(3) x Z{2) formu- 
lation with periodic boundary conditions. The dynamical variables are 
constrained by the required coincidence of Z(2) and SO (3) monopoles. 
We propose an update algorithm that satisfies the constraints and is 
straightforward to implement. We further prove that this it reaches all 
configurations. We show how the boundary conditions put constraints 
on the configuration space. We measure gauge invariant vortex coun- 
ters for "thin" , "thick" and "hybrid" vortex sheets. For comparison we 
also measure projection vortex counters defined in the maximal center 
gauge. 



PACS indices: 11.15.Ha, 11.30.Ly. 



1 Introduction 



Lattice QCD continues to maintain an important role in the search for the physics of color 
confinement. The lattice regulator maintains gauge invariance at all costs. Dynamical 
variables are group elements rather than elements of a Lie algebra. As a consequence 
many of the topological features that are prominent candidates for elucidating the physics 
of confinement have natural lattice definitions. These include U(l), Z(N), SU(N)/Z(N) 
monopole loops, Dirac sheets, Z(N) and SU(N) vortex sheets etc. These objects are often 
abundant in U(l) and SU(N) lattice gauge theories. They become singular only as one 
approaches the continuum limit. 

Consider the case of SU(N). Further consider a multiply connected region in which 
all links are gauge equivalent to 1 on any simply connected patch of the region, i.e. all 
neighboring plaquettes = 1. Then the value of a Wilson loop lying in this region would take 
the value of a center element of SU(N) = e 2mn ' N , which for indicates the trapping of 

a vortex. The occurance and absence of vortices gives a fluctuating value that can disorder 
the Wilson loop and lead to an area law. 

Yaffe|jl|, Tomboulis@ and Kovacs and Tomboulis0 have developed a formulation of 
SU(N) gauge theory that is manifestly SU(N)/Z(N) invariant. In this formulation, center 
elements, Z(N), multiplying each link leave the action and measure invariant. New Z(N) 
variables, defined on plaquettes, a(p) , carry the Z(N) degrees of freedom. This formalism, 
equivalent to the standard SU(N) form, allows an elegant topological classification of the 
SU(N)/Z(N) and Z(N) vortex configurations occuring on the lattice. See also references 
[4-7]. 

In this paper, we address the issues of simulating the Tomboulis formulation on a 
periodic lattice. Many results have followed from this formulation without doing simulations 
in these variables^]. As a first calculation, we tag Wilson loop measurements by the 
occurance of vortices linking the loop using a number of different vortex counters. We 
restrict our attention to SU(2) in this paper. 

We also measure the P (projection) vortex counter in the original SU(2) formulation 
following [9-17] for comparison. Projection vortices arise in a Z(2) gauge theory derived 
from the original SU(2) theory by going to the maximal center gauge and then replacing 
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links by sgn(tr(L r XiAt )). The projected theory has "thin" Z(2) vortices defined on one lattice 
spacing. They have been found to be well correlated with center vortices and therefore a 
measurement of P vortices is a predictor of them. 

The thrust of this paper is to study the configuration space in the Tomboulis variables 
on the torus. These variables are subject to constraints in a rather indirect way. We propose 
a constructive update algorithm which is straightforward to implement and we prove that 
it reaches all configurations, i.e. it is ergodic. 

In section 2 we rederive the Tomboulis form of the partition function. There are config- 
urations on the torus that give zero weight which we exhibit in section 3. In section 4 we 
elucidate two alternative definitions of the configuration space for the Z(2) variables cr(p), 
the indirect definition and constructive definition. Appendix B gives the proof of their 
equivalence. In section 4, we discuss various vortex counters. In section 5, we measure 
these vortex counters for Wilson loops. This allows us to tag Wilson loops to study the 
disordering mechanism. In Appendix C we consider the case of anti-periodic boundary 
conditions [0. 



2 Derivation of Z Z (2)xSU{2)/z(2) on a torus 

Our derivation here is due to Tomboulis @, and reviewed recently by Kovacs and Tomboulis 
[[| for free boundary conditions. The latter paper gives a thorough pedagogical review of 
the formulation and the topological features. In addition they have paid close attention to 
visualization of vortices, both as surfaces and as curves in 3-d slices. 

Consider the partition function for the Wilson action, 

Z = J [dU{b)\e* V (^Y,\^U{d P )]}. (1) 

Define a Z(2) variable on the links (bonds), 7(6), and insert a constant into the partition 
function, 

Z = £ / W{b)l{b)} exp U £ hr[U(dp)]) • (2) 

7(6) \ p / 

(See Appendix A for Z(2) notation, algebra, characters, delta functions, etc.) Apply the 
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Haar invariant transformation 



U(b) -> U(bh(b); Z = Y.J l dU ( b )} ex P te^W)W*)) • 

7(6) V p / 

Isolate the sign of the plaquette, rj(p), tv[U(dp)] = \ti[U(dp)]\ x rj(p), writing 

z = E / ex P (/?E ^MtfWlMSptofr)) • 
7(6) \ p / 

Next introduce a new Z(2) variable, <r(p), defined on plaquettes, 

1 = ^i(a(p) x -y(dp)r)(p)), 

<r(p) 

Z = W W)] E%(P)X7(^(P)) exp f)9E^N^)]k(p)V 

7(6) Uw J v p / 

Expand the delta function in Z(2) characters 



S(a(p)rj(p) x 7 (9p)) = - E (xt( p ) (cr(p)v(p)) Xr(p) 

T(p) 

= \ E {xr( P ) (<r(p)77(p)) x 7 (9p) • 
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T(p) 

This gives 

7(6) <t(p) r(p) 

{LI \xrip) HpMp))} |nX7(0 P ) (r(p))| exp ^\HU(dp)}\a(p^ . 

In the second character, we can rearrange the product over plaquettes into a product over 
links, 

lix^dp) (r(p)) = 11X7(6) ( r (<*)) • 

P b 

Now do the 7(6) summation 

z = /WOIE 

<r(p) 

En^w K^w)n^(^ 



2 

r(p) P - 



<r(p) 



expf/^E^NC/^)]!^)]. (3) 



p 



The item in square brackets is the starting point for much of the analysis in this paper. 
The dynamical variables are {U(b),a(p)}. The r(p) sum can be done. 



C(a(p)r](p)) 



EllXrW (a(p)r](p))l[6{T(db) 

r(p) V b 



(4) 



Y[ 5 {a(dc)rj(dc)) x <{ . ^ x constant (5) 

The constraint 5 (r(db)) means that there must be an even number of r = — 1 plaquettes 
in the co-boundary of the links (i.e. the six plaquettes contiguous with the link). 

The last equality needs further explanation. If a cube contains an odd number of faces 
with 7] = — 1 then by definition it contains an SO(3) monopole. Similarly if a cube contains 
an odd number of faces with a = — 1 then it contains a Z(2) monopole. The delta function 
on the cube requires that any SO (3) monopole be paired with a Z(2) monopole at the 
same location. We show in Sec. 4 that for the integrand of the partition function to be 
different from zero it is necessary to have Y\ C S (a(dc)r](dc)) = 1. However on the torus, 
this is not sufficient. There are configurations, {U(b),a(p)}, for which the delta functions 
on the cube are unity yet the integral vanishes. We denote these configurations " weight = 
0." In section || we give examples of such configurations which are co-closed vortex sheets 
that wrap around periodic boundary conditions. We further show that when the r sum 
differs from zero, it is a constant in the variable 

a(p) = a{p)r]{p). (6) 

We denote these "weight = 1" configurations. In section ^ we exhibit an update algorithm 
that reaches all weight = 1 configurations and respects all constraints. 

Restricting {U(b)} and {a(p)}, to the weight = 1 configurations we obtain 

Z= £/ [dU(b)]'l[8(a(dc) V (dc))exp U^\\tr\U(dp)]W(p)) ■ (7) 

<r(p) c \ P / 

Note that this form is invariant under U(b) — > ^(b)U(b). All configurations related by this 
transformation are SU(2) representatives of the invariance group SO(3). The Z(2) part is 
explicit in the a(p) variables. 
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3 Zero weight configurations on the torus 



In this section we find configurations on the torus which have zero weight as indicated in 
Eqn.fEp. We take periodic boundary conditions in all directions. The zero weight reflects 



the fact that a vortex wrapped around the torus is topologically stable p8fl. It can not be 
reached from a non-vortex configuration. This property is inherent in the formalism. 

Our starting point is Eqn.(^). Using Eqn.(^) this becomes 

C(a(p)) = EllXrW (a(p))i[5 (r(db) • (8) 

r(p) V b 

The t delta function constraint requires that the plaquettes forming the co-boundary of 
any link must occur in even numbers. Clearly a closed surface made by tiling with r = — 1 
plaquettes, with r = +1 elsewhere, will satisfy all these constraints. 

Now lets turn to the a variables. Consider a configuration in which a(p) = — 1 on all 
plaquettes pH ke and = +1 elsewhere. The upper indices denote the plaquette orientation, 
and the lower indices are the space-time coordinates, (i,j,k,£). This configuration is a 
co-closed vortex sheet wrapped around the torus in the 3 and 4 directions. If a(p) = =pl 
and r](p) = ±1 everywhere on this co-closed vortex, it is a cr/i] vortex. If these two cases 
occur on different patches of the vortex sheet, then it is a hybrid vortex with co-closed 
monopole loops at the boundaries between the patches. 

The a cube delta function constraints, Eqn.(^|), are satisfied because cubes will either 
have no vortex plaquettes in common, or will have two vortex plaquettes on opposite faces. 
In spite of this we now show that this a configuration has zero weight in the partition 
function. 

Consider Eqn.(|8]) applied to an arbitrary function of F(r(p),a(p)). Define the set C 

t(jj) b C 

i.e. it is the set of all r configurations with the property that l\ b 5 (r(db)) = 1. 
The set C forms a group under the multiplication defined through 

(r 1 T 2 )(p) = r 1 (p)r 2 (p). (9) 
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Given that T\ and r 2 are group elements then each has an even number of negative plaquettes 
in the co-boundary of any link. Clearly the product will have the same property. The 
identity element is r(p) — 1 for all p, and the elements are their own inverses. 

Using the invariance property of the group summation: 

^F(r r;a) = ^F(r;a), 

where r is any element of the group C. 
Substituting for F our case reads: 

C(«(P)) = J2IiXr ip) (a(p))l[6(r(db)), (10) 

r(p) P b 

= E II Xr (p)r(p) (a(p)) I] 5 ( T (^>)) > 
r(p) P 6 

= nX7 ( P )(a(p)) x C(o!(p)), (11) 
p 

for any r G C Therefore if we can find a group element r for which Y[ p Xt ( p ) ( a ip)) 1 
then C(a) = 0. 

Choose r (p) = — 1 on the subset of plaquettes: {p}f 00 } for all i,j and +1 elsewhere. 
It is obviously a member of the set C. This is a closed tiled surface of 1,2 plaquettes 
wrapping around the 1, 2 directions for the 3 and 4 coordinates fixed to 0. Using the fact 
that — X+(— ) — X-(+) = +1 = — X-( — )> we count the number of sites where both 

a — — 1 and r = —1. 

The configuration under consideration is a co-closed set of 1, 2 plaquettes wrapped 
around the 3,4 directions and the r{p) configuration is a closed tiled set of 1,2 plaquettes 
wrapped around the 1, 2 directions. They have only one negative plaquette in common at 
position (0,0,0,0). Therefore I\ p Xt (p) ( a (p)) — — 1 implying C(a) = 0. 

4 a(p) configuration space 

We have seen examples in the last section of the interplay between the r and a configu- 
rations in finding non-zero contributions to the partition function. We are interested in 
simulating in the variables {U(b), &{p)}. Thus far the allowable a (= cr x rj) configurations 
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are defined indirectly in terms in the allowable r configurations which are also indirectly de- 
fined by constraints. The corresponding simulation would also be very indirect and perhaps 
difficult to implement. 

We propose a constructive definition of allowable a(p) configurations by building them 
up from "star transformations", i.e. correlated sign flips of the a plaquettes occurring in 
the co-boundary of each link. Since these are constrained updates of six plaquettes, it is 
not clear that we can reach all allowed {a} (or equivalently {a}) configurations by this 
method. However we show that this definition of allowed a(p) configuration is identical 
to the above definition. The proof is relegated to Appendix B. In this section we give a 
summary of the result. 

Before discussing the a configurations let us first describe the link updates. This is a 
straightforward generalization of the link updates for SU(2). The proposed change in a 
link might change the sign of the rj plaquettes in the co-boundary of the link. If one of 
these changes sign, we need to flip the sign of the corresponding a plaquette so that the a 
configuration is unchanged. Then the Monte Carlo step is essentially the same as for the 
SU(2) update. 

Next consider the above mentioned "star transformations". Our proposed update is to 
flip the sign of the six a plaquettes forming the co-boundary of the links. 

Assume we are starting from a weight = 1 configuration, Eqn.([5]). It is easy to see that 
both these update steps will preserve the cube constraints. 

In the previous section we described vortex sheets that wrap around the torus. Consider 
the operator constructed out of /i, v plaquettes || 

= II v(p)°(p) = viS^viS^) = ±1, 

where S^ u is a whole tiled fi, v plane. M^ v = ± for an even/odd number of vortices of 
stacked fi, v plaquettes wrapping around the orthogonal directions of the torus. 

We start the algorithm with this M^ v = +1 in all 6 planes. It is easy to see that our 
update algorithm preserves N^ v . Hence, starting with non-zero weight configurations, we 
do not generate the zero weight configurations described in the last section, since this would 
involve AAf^ tU ^ 0. 
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Let us define the relevant sets of configurations more carefully. Consider Eqn. fllOD 
CV) = EIIXr W (a(p))n*(^)). 

r(p) P 6 

= X)IIXT(p)(a(p)), 

The third line is a shorthand for the product over characters (see Appendix B). The 
bracket, (a, r) = ±1, is negative if and only if there are an odd number of plaquettes for 
which both a(p) = —1 and r(p) = —1. 

The second line is an alternative way to specify the sum, where: 



C = {reA\l[S{r(db)) = l}. 



Configurations form a group, Eqn.([|). C is a subgroup of the group A of all configura- 
tions, |*4| = 2 6N in number, where N is the number of lattice sites. Restating: 

C is the group of all configurations {r} with an even number of r — — 1 pla- 
quettes occurring in the co-boundary of every link, i.e. forming a closed tiled 
surface of negative plaquettes. 

There is a second group of interest, 

C = {a E A\(a,r) = 1 Vr 6 C}. 

Restating: 

C is the group of all configurations {a} for which C(a) may be different from 
zero. We will show that on this set we have indeed C(a) ^ 0. (See appendix 
B.2, Proposition 2) Recall from Eqn. (j^) that if we can find a single configuration 
t for which (a, r ) ^ 1 then C(a) = 0. 

Therefore C is the group of a configurations which have weight = 1 in the sense of 
Eqn. ([5D, i.e. the configurations that contribute to the partition function. Further, C is 
the group of r configurations that form closed tiled surfaces as required by the explicit 
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constraints in Eqn. ([l2|) . In the previous section, we found a zero weight configuration a by 
finding a configuration r for which (a,r ) = — 1. 

The group C has only an implicit definition here. The group C has an implicit definition 
in terms of this C. Therefore its definition is even more indirect. Even without an explicit 
definition, we have been able to specify precisely those configurations {a} that contribute 
non-zero weight to the partition function. 

There is a third group of interest, 

V = { a eA\a = ]Ja Bbt }. 

where refers to an individual star transformation on the z'th link bi, and the product 
indicates all possible products of them. 

Restating: 

T> is the group of all configurations {a} which can be built out of products of 
"star transformations" starting from the identity configuration. 

This is the constructive definition that is straightforward to implement in a simulation. 
The proof in Appendix B shows that the group T> is identical to the group C. In this way we 
have shown that by our proposed algorithm is ergodic, reaching all configurations allowed 
in Eqn.(||). 

Let us return to Eqn.(|) and the "cube constraints." Using our definition, a(p) = 
g(p)t](p), the constraints are written as Y\ c S(a(dc)) = 1. The cube constraint simply 
asserts that for a configuration to give a non-vanishing contribution, every cube in the 
lattice must have an even number of faces with a = — 1. 

Let us suppose that a particular cube has an odd number of faces with a = — 1. Then 
consider a configuration r which takes values —1 on all 6 faces of this particular cube. 
This is a closed surface and therefore satisfies the constraints imposed on r. For this case 

HXroip) («(?)) = -1 
v 

and therefore by Eqn. (|TTD C(a) = 0. 

The reason for signaling out this necessary constraint is that it is local. The zero weight 
configurations described in the last section necessarily wind around the torus. 
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5 Simulation of vortex counters 



Vortices have long been considered as prime candidates for the essential dynamical vari- 
able to describe confinement. A simulation offers a tool that allows one to correlate the 
occurance of vortices with values of other dynamical variables. Hence as a first application 
we use this formalism to measure various vortex counters for Wilson loops. 

Consider the SU(2) formalism with the standard Wilson action. Further consider a mul- 
tiply connected region in which all links are gauge equivalent to 1 on any simply connected 
patch of the region, i.e. all plaquettes = 1. This could occur if the vortices are very dilute 
and have a cross section small compared to average separation. 

Then the value of a Wilson loop lying in this region = ±1, the center of SU(2), corre- 
sponding to an even/odd number of SU(2) vortices linking the region. The occurance and 
absence of vortices gives a fluctuating value that can disorder the Wilson loop and lead to 
an area law. 

In the £0(3) x Z(2) formulation the Wilson loop is given by 

W[C] = (^tr[C] Vs a s } c =8S, 

where r/s and as are products of rj and a over any spanning surface Jj], |], j3j. 

Kovacs and Tomboulis|3[] define three vortex counters for thick vortex sheets, thin vortex 
sheets, and hybrid (patches of each on the sheet). Their definitions require measurements 
on all spanning surfaces. We measure here only the minimum spanning surface. Hence we 
must interpret our measurements as best we can in this limited simulation. 

• Thin: 

•A/thin = a s- 

If this value, = ±1, is independent of the spanning surface, then this counts thin 
vortices. 

• Thick: 

^thick = Vs sga[trW(C)}. 
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This object is counting something more elusive since unlike the above case, the vortex 
structure is spread over many lattice spacings. Nevertheless it is always possible to 
find a representative of SO(3) such that the r\ vortex defines the topological linkage^, 
18|j . The f] vortices can be deformed by a Z(2) transformation of links giving different 



representatives of SO (3) without cost of action. One can move a linked r] vortex sheet 
so that it no longer links the Wilson loop and further even transform it away However 
in this case the negative contribution will be transferred to one of the perimeter 
links of the Wilson loop, and it will not affect the value of A/t^ick- Again if this is 
independent of the spanning surface, then this counts "thick" vortices. 

• Hybrid: 

^hybrid = ^thin x ^thick = <*SVS sgn[trW(C)]. 

As one considers all spanning surfaces, the sign of A^hm might change. However if 
the sign of A/thick arwa y s compensates then this counts hybrid vortices. 

We measure these three counters for Wilson loops, taking the minimal spanning sur- 
face. Then, for example, A/thin = ~~ 1 does not distinguish thin from hybrid, and similarly 
A/"thick = — 1 does not distinguish thick from hybrid. However they do measure the oc- 
curance of an object piercing the spanning surface responsible for sign fluctuations of the 
Wilson loop. Notice that A/Hybrid is just the sign of the Wilson loop, Eqn. fll3|) , and can be 
calculated in the original SU(2) formulation. 

A fourth definition of vortex counter enters in the work of Del Debbio, Faber, Greensite 
and 01ejnik[|. See also Refs.(|L0|, 0, 0, 0, 0, [15|) Starting with the SU(2) formalism, 



go to the maximal center gauge which maximizes 

Then consider a Z(2) gauge theory in which the links are replaced by their Z(2) values 

sgn(tr(f4 )At )). 

Denote the plaquettes in this Z(2) gauge theory by £(p). The negative plaquettes of this 
theory, = —1 form thin vortices, residing on one lattice spacing similar to a vortices. 
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They are denoted P (projection) vortices. The authors find evidence for a strong correla- 
tion between P vortices and thick objects, denoted center vortices, analogous to the thick 
vortices of Tomboulis. Their calculations proceed in the original SU(2) theory with the 
added observable of the P vortex counter to segregate contributions to Wilson loops. 

• Projection: 

•^projection = £s- 
This object is independent of the spanning surface. 

We also measure the P vortex counter for comparison, using the SU(2) formalism. 



6 Numerical Results 

Simulations were done on a 12 4 lattice for (3 = 2.30. Measurements were binned to 10 bins 
and jackknife errors calculated, [thin]: 200 ; [thick]: 400 ; [hybrid]: 200 ; [projection]: 1000 
measurements. We monitor the coincidence of Z(2) and SO (3) monopoles which can slip 
due to round off error. 

Let us emphasize again that in our application here these counters measure only the 
occurance of various objects piercing the minimal spanning surface, not the species of 
vortex. More general analysis will be given elsewhere [T3 . 



Fig. 1 shows the fraction of Wilson loops which have an even number of vortices, x e , 
as a function of the area. All the counters approach 50% from above with similar behavior 
(x = 1 — x e ). They are each counting different things and are not expected to be equal. 

A/hybrid is a good reference curve since it measures the sign of the Wilson loop itself. 
The area law arises from a near cancelation of fluctuating values due to approximately 
equal occurance and absence of thin, thick or hybrid vortices. 

The "thick" curve lies on top of the "hybrid" one. Hence the added factor of a,s in 
the hybrid counter has little effect here. We expect a vortices to be heavily suppressed for 
increasing f3 since they cost action proportional to the vortex area. However aX (3 = 2.30 
the density of Z(2) (or SO(3)) monopoles = 0.2155(2) (Random plaquette signs would give 
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5 10 15 20 25 

Area 

Figure 1: Fraction of Wilson loops with an even number of vortices piercing the minimal 

surface, x e . (x = 1 — x e ) 

a density of 0.5). Hence in spite of the near coincidence of these two curves, a vortices 
and a patches of hybrid vortices are important at this value of f3. There is further evidence 
below. 

Since Wilson loops are positive, one expects that the fraction of even loops x e should 
always dominate, as they do in all cases. 

There are two special points on these curves. A/thi c k = +1 by definition for a plaquette, 
i.e. a single "thick" vortex can not link a 1 x 1 Wilson loop. Hence the fraction x e = 1. 
Further this gives the fraction x e equal for the two cases A/thin, and A/Hybrid 

The AAprojection case gives the same result as reported in Ref 0. 

Figs 2-6 shows Wilson loops, W, and the tagged Wilson loops corresponding to even: 
W e and odd: W number of vortices as a function of loop area and their logarithmic 
derivatives. 

• Fig. 2: [Thick or thin segments piercing the minimal surface] This vortex counter 
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Figure 2: Wilson loop, W, and tagged Wilson loop with even or odd number of "thin" 
or "thick" vortices piercing the minimal surface. 

is just the sign of the Wilson loop itself. The positive and negative contributions to 
the Wilson loop are averaged separately. Since each component has approximately 
an equal and opposite asymptote, this illustrates the cancellations due to disordering 
and the difficulty of measuring large loops. 

• Fig 3: [Thick segment piercing the minimal surface] By definition W e has an even 
number of thick segments piercing the minimal surface. Yet it still has an exponential 
fall off with area. See Fig. 6 which gives the logarithmic derivative showing that W e 
has about half the string tension of the W. The thin segments are still present and 
they account for the disordering. 

• Fig 4: [Thin segment piercing the minimal surface] By definition W e has an even 
number of thin segments piercing the minimal surface. It nearly coincides with W 
and hence cancellations due to thin segments not important and are lost in the noise. 
This curve also shows a breakdown of the correlation of the sign of W a and the sign of 
the vortex counter for all but the small loop sizes. This object is not a good predictor 
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Figure 3: Same as in Figure 2 but tagged by "thick" only, 
of the sign of the Wilson loop. 

Fig. 5: [Projection vortices piercing the minimal surface] For loops of area 9 and 
higher, the sign of the vortex counter correlates with W a . The last four points coincide 
with those reported in Ref. 0. Comparing Fig. 3 [thick] and Fig. 5b. [projection] the 
latter data are about a factor of 10 smaller. 

Fig. 6: [Logarithmic derivatives] W and W e [thick] show a constant string tension for 
larger loops. We suspect that with better statistics, W e [thin] will also. W e [hybrid] 
shows the vanishing of the string tension if one removes the disordering mechanism 
completely. Larger loop areas are needed to decide if the logarithmic derivative of W e 
[projection] will also go to zero, or stablize which would indicate that a disordering 
mechanism remains. 
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Figure 4: Same as in Figure 2 but tagged by "thin" only. 
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Figure 5: Same as in Figure 2 but tagged by "projection" only. 
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Figure 6: String tension: K = — ^ d ^ ea , similarly for K e . 

7 Conclusion 



We have proposed a simulation algorithm for the partition function in the Z(2) x 5*0(3) 
formulation. We show that the algorithm is ergodic, reaching all relevant configurations. 
We consider vortices which wrap around the torus. We find that these have zero weight in 
the partition function reflecting the fact that they are topologically stable. 

Thick vortices are known to be an important factor in disordering the Wilson loop. 
As a first calculation, we measure various vortex counters, in order to see how they are 
correlated with other observables. 
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A Appendix on Z(2) Algebra 



We follow the definitions of Tomboulis and Kovacs [0, [| We used b, p and c to denote the 
links, plaquettes and cubes respectively. Occasionally we use indices (e.g. v%ki) t° denote 
the objects location (ijkl) and orientation The d and d operators have the usual 

meaning: the boundary and coboundary operators. 

Z(2) is the multiplicative group with two elements and we will denote it's elements 
with Greek letters: a, (3 ... The Z(2) group admits two representations: M + (±l) = +1 
and M_(±l) = ±1. The characters of the representations are: = +1 for M + and 

X-i(±l) = ±1 for M_. 

The Z(2) delta function is defined as: = 1 and 5(— 1) = 0. 

We list the basic of the properties of the characters: 

Xr(a(3) = Xt(u)xAP), 
6 ( a ) = 2^ Xr ^- 



B Appendix, Proof of ergodicity of the a update al- 
gorithm 

In this appendix we prove that the two groups described in Sec. 4: T> built up by "star 
transformations" on each link is identical to the group C defined by constraints on each link. 
Sections B.l, B.2 and B.3 give preliminaries. Section B.4 proves the result by induction. 
We define intermediate groups T>e built out of star transformations on a subset of links 
and similarly Ce restricted by constraints on the same subset of links. We then increase 
the set to E 1 by an additional link and proceed by induction. 

The partition function contains the factor, Eqns. (0, §, §, |, 0, [12|) 

=EIIXT(p)(Qt(p)) = J! (a, r). 

We will see that C[a] = for a ^ C and C[a] = \C\ = number of elements in the set C for 
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a G C, where 

C = {a G C\(a,r) = 1, Vr G C}, 

and where C is defined to be: 

C = {r6 *4| n^( r (^)) = !} ^ C = {r G *4|r(<%) = 1, V6 G B}. 
b 

and where S is the set of all links of the lattice. 

Recall C is a closed tiled surface of negative r plaquettes as required by the constraints 
in the partition function. The r variables are summed, leaving the a variables. C contains, 
e.g., a vortex of stacked negative a plaquettes which have non- vanishing weight in the 
partition function. 

B.l Notation 

P denotes the set of all plaquettes of the lattice and B the set of all links. A "configuration" 
is defined to be a function that associates an element of Z 2 to each plaquette. We will denote 
the configurations with Greek letters and write: 

a : P — » Z 2 p i— > a(p) G Z 2 . 

The set of all configurations is denoted by A = {a : P — > Z 2 }. For every set K we will 
denote with \K\ the number of elements in the set. \A\ = 2 6N where iV is the number of 
sites in the lattice. 

Configurations form a group under multiplication: let a and (3 be two configurations in 
A. Then a[3 is defined to be: 

a f3{p) = a(p)(3(p). 

We denote the unit element 1 G A which assigns +1 to all plaquettes. Elements are their 
own inverses. 

The partition function, Eqn. |3|, contains a summation over all r configurations that 
satisfy the constraints: l\ b 5 (r(dbfj = 1. We denote this set by C: 

C[a] = EIIXr(p)(a(p)), (13) 

rec p 
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where C = {r e A\ U b $ (^(db)) = 1}. 

Simplifying the notation in the summand we define: 

(r,a) = Hxt(p)(u(p)). 
v 

We will list here some properties of the bracket (•, •) without proof. Lets take a,/?, 7 € A. 
Then: 

(a, (3) = ((3, a), 
(aft, 7) = (a,7)(/3,7), 
(a,/37) = (a,/3)(a,7), 

(a,l) = (l,a) = l. 

B.2 Some theorems regarding the subgroups of .4 

Consider an arbitrary subgroup of ^4 denoted /C, (of which C is an example). 
Proposition 1: Let /C be a subgroup of A and define: 



Then: 



K[a} = (a,p )K[a] (14) 



for any f3 e /C. 

Proof: Since /C is a group we have: 

£ (a, /3 /3) = £ <«>0> 
/3e/c /3e/c 

where we used the property of the group sum that {PoP} is a rearrangement of the group 
elements Using the properties of the bracket we have: 

K[a] = £<a,/?o/?> = £<a,A,><a,/?>, 

= (a, fa) («, P) = (a, A))-K"H- 
/3e/c 
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Definition: Let /C be a subgroup of A. We define: 

K = {a E A\{a,P) = 1 V/3 e /C}. 

£ (e.g. C) has always at least one element, 1, since (a, 1) = 1 for any a E A. Moreover it 
is easy to prove that K is a group too. Now we note the following lemma: 

Lemma 1: A — {1}. 

Proof: Choose a E A with a ^ 1. This means that there is at least one plaquette 
Po E P for which a(p ) = —1. Then if we take (3(p) = 1 for all p ^ p and f3(p ) = — 1 
which is an element of A we see that a and (3 have only one plaquette, p , on which both 
are —1. Then 

(a, (3) = YlXa(p)(P(p)) = -1, 
p 

proving that a £ A. Thus we proved that if a ^ 1 then a ^ A Hence .4 = {1}. 

Proposition 2: Let /C be a subgroup of A. If if [a] is defined as in proposition 1 we 
have K[a] = for a E~ IC and K[a] = \K\ for a E IC. 

Proof: Using proposition 1 we have: 

K[a] = (a,(3 )K[a\ 

for any /3 E IC. Clearly if (a,/3o) ^ 1 then K[a] = 0. Thus K[a] = for all a that have 
at least one element (3 E JC for which (a,f3 ) ^ 1. This is equivalent with saying that if 
a £ K, then K[a] = 0, thus proving the first part. 

Now consider a E IC. This means that (a, (3) = 1 for all (3 E IC. Then: 

This concludes the proof. 

Theorem 1: Let IC be a subgroup of A. Then |/C||£| = |^4|. 
Proof: Consider the following: 

I=J2 K[a] = \K\ 1 = l^ll^l, 
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where we used proposition 2. Interchange the sums in I and using the commutative property 
of the bracket: 



i = E K\a\ = £ 5>,/?> = ££</?,«> = E m- 

where A\J3] is defined exactly as K[a] using the fact that the bracket is commutative. Now 
proposition 2 tells us that A[a] = \A\ for a E A and otherwise. Using lemma 1 we have 
that A[l] = \A\ and zero otherwise and thus we get I — \A\. But since we already know 
that / = |/C||^| we have: 

\£\\ic\ = \A\. 

B.3 Constructing V 

We already know that C is a subgroup of A. In this section we construct a subgroup of C, 
denoted V, which, in the next section, will be shown to be the entire C. 

Definition: We will call "a star configuration around link b" the following configura- 
tion: 

' +1 p £ db, 
. -1 pe db. 



Lemma 2: For any link b we have aQ b G C. 

Proof: We can prove that ag b is a member of C by checking that: 

(a db ,r) = l, WreC. 
Consider any link b and any r G C and compute: 

( a db' T ) = I[Xa Sb (p)(r(p)) = [] Xi(r(p)) I] X-i(r(p)) = X-i(r0b)). 

p pgdb pedb 

Since r G C then r(db) = 1 for any link b. Then: 

(a 6b ,T)=x-Mdb))= X -i{l) = +l. 
for all r G C. This proves that for any link b we have ag b G C. 
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If we take the set of all star configurations by taking all possible products between them 
we can generate a group, T>. Since all the star configurations are included in C, which is a 
group itself, the group V is a subgroup of C. Lets write V explicitely: 

V = {aeA\a = l[a Bbi }. 

where the product is over all configurations reached by star transformations. In order to 
prove that T> covers all C we will show that T> has the same number of elements as C. Since 
V C C we have \V\ < \C\. All we have to prove now is that \V\ > \C\. 

B.4 Proof that C = V 

We now prove that \V\ > \C\. Using theorem 1 we know that = |»4|. Thus we can 

prove that \T>\ > \C\ by proving that \C\\V\ > \A\. Lets recall the definitions of these sets: 

C = {re A\r0b) = 1, \/b G B}, 
V = { a e A\a = Ua Bbi , b t e B}. 

As we see the set C is constructed by means of eliminating the configurations that do 
not obey a certain constraint whereas the set V is constructed by generating all possible 
combinations built from star configurations acting on the identity configuration. These sets 
admit a generalization as follows. Define: 

C E = {re A\r(db) = 1, Wbe E}, 
V E = {aeA\a = '[[ag bi , h e E}, 

where E C B. It is easy to check that both C E and T> E are groups. In words, C E is the set 
of all the configurations that obey the constraint only on the subset E of all links and V E 
is the group generated by star configurations associated only with links in E. It is obvious 
that C B = C and that V B = V. If we prove that |X>e||Ce| > |»4| we will implicitly prove 
that \V\\C\ > \A\ and thus proving that \T>\ > \C\. 

Theorem 2: \V E \\C E \ > \A\. 

Proof: We will prove this using induction. The first step will be to prove this for E = 
and the second step will be to prove that the relation holds for E' = E U {&}, where b is 
any link, assuming that the inequality holds for E. 
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B.4.1 Initializing 

Let E — 0. We have = A since there is no constraint that has to be obeyed and Vq, = {1} 
since this is the group that has no star configuration in it. Then we have: 

|P ||C | = lx |^| = |^|. 

Therefore the inequality is verified for E — 0. 

B.4.2 Iterative inductive step 

Assume that: 

\T> E \\C E \ > \A\ 

for a certain E. Let us prove that this also holds for E' — E U {b} for any b ^ E. We will 
see what happens with Ce and V E when we increase E by one element. 

Iteration on Ce 

When we increase the number of elements in E, Ce grows smaller since there will be 
more constraints to obey. It may happen that the new constraint, namely r(db) — 1, is 
superfluous i.e. all configuration in Ce already satisfy this constraint. In this case Ce> = Ce 
and \Ce'\ = \Ce\- Lets now see what happens if there is at least one element in Ce that 
doesn't obey the new constraint. In this case we can break down Ce in two disjoint sets: 
Ce' the sets of all configuration r e Ce that obey the new constraint and R the set of all 
configurations r e C E that do not obey the constraint. In case that R lets take r e R. 
Using this we can construct a one to one mapping between Ce> and R: 

f:C E >^ R, f(r) = t t, 
f-^.R^Cv, f(r)=r r. 

It is easy to check that / is indeed one to one. 

Now since there is a one to one mapping between Ce 1 and R we have that \Ce<\ = \R\- 
But since they are disjoint sets of C E and C E ' U R = C E we have \C E >\ + \R\ = \C E \- Thus 
we have \Ce>\ = \\Ce\- 
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Summing up we know that by adding a new constraint \Ce>\ is either equal with \Ce\ 
when the constraint is superfluous or is \\Ce\ when we have at least one element in Ce that 
violates the constraint. 

Iteration on V E 

Now we will look at V E '- What happens when we pass from V E to V E >7 We add a 
new element in the group. Thus the group grows larger with one exception: it may happen 
that ag b G T>e although b G" E. Then the group will stay the same since every combination 
that doesn't involve ag b is already in the group V E and every combination that involves 
a db can be generated by elements already in the group. Formally, if a db G T>e then for 
any a G V E we have a db a G V E and thus V E > = T> E - This proves that if a db G T>e we 
have \T>e>\ = [De\- It is not obvious that this case happens but we can easily construct 
one: take a site in the lattice and take the eight links that form its coboundary. Since the 
coboundary of a coboundary is nul (i.e. a 2 = 0) it means that if we perform all the star 
transformations associated with the links in the coboundary of the site we get the identity. 
This means that the product of any seven of this star transformations is equal to the eighth 
star transformation. 

Consider what happens when a§ b ^T>e- Then for every element a G T>e we have two 
elements in V E >'- a E V E > and ag b a G T>e'- We can actually construct two sets in V E '- 
T> E which is trivially included in V E ' and ag b x V E = {« G T> E >\oi = ®Q b P, P G T> E }- If 
a db G" T> e then these two sets do not overlap. 

If they were to overlap then there is a G V E H ag b x V E with a G V E and a = ot db f3 
with p G T> E - Then a db = aP and since both a and P are members of V E then ag b G V E 
which contradicts our assumption. 

It is easy to see that T>e and a:^ x T>e have the same number of elements and that 
T>e U ag b x T>e = T>e>. Now since T>e U a db x T>e = T>e' and T>e fl ag b x T> E = we have 
\D E \ + \ag b x T> E \ = \V E '\. Moreover \V E \ = \ag b x V E \ and thus \V E '\ = C 2\D E \- 

Summing up if ag b G T>e then \T>e'\ = [De\ and if G" we have = 1\De\- 

Intermediate summary of possible cases 

Thus we arrived at the conclusion that if we add another link b to E one of the four 
following things may happen: 
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• the new constraint is superfluous (\Ce'\ = \Ce\) and a 9b G" V E (\T^E'\ = 2\Ve\). Then 
\C E '\\V E >\ = 2\C E \\V E \ > \A\. 

• the new constraint is superfluous (\C E >\ = \Ce\) and a 9b G V E (\T>E'\ = \^e\)- Then 
\Ce>\\T>e>\ = \C E \\T> E \ > \A\. 

• the new constraint is not superfluous (\Ce>\ = §|Ce|) and a 9b G" T>e {\T^e'\ = 2|X>e|). 
Then \C E '\\V E ,\ = \C E \\V E \ > \A\. 

• the new constraint is not superfluous (\Ce>\ = and a 9b G V E (\T^E'\ = I^bD- 
Then |Ce/||De/| = ||Cb||X>e] and we don't know if this is smaller or greater than \A\. 

If we look at the scheme above we see that if we can prove that the last case never happens 
than we proved our theorem. This is exactly what we prove in the following lemma: 

Lemma 3: If a 9b G V E then the new constraint is superfluous. 

Proof: If a db G V E then there is a subset E C E with the property ag b — Hb'eE a §b> 
(this simply asserts that a db can be written as a product of star configurations associated 
with a subset of links in E). 

Lets denote with P the set of all plaquettes that form the coboundary of E : P = dE 
All the plaquettes in P have at least one neighboring link in E and thus they are flipped 
at least once when you take the product Tib'eEo a db>- 

Since the whole product is equal to a db the coboundary of b, db, has to be included in 
Pq. Otherwise we cannot flip the plaquettes around b and we cannot form a db . 

Further since the final state is a 9b we know that all plaquettes p E P — db are flipped 
an even number of times and all plaquettes p G db are flipped an odd number of times. 

We define a function k(p) on plaquettes in P Q that returns the number of times this 
plaquette is flipped. Using this it is easy to prove that n{p) is even for p G Pq — db and 
odd for p G db: 

Yl a db ,(p) = J] = (-l) E "'^o e »' (p) = (_1)"(p) 

b'eE b'eE 

where ey{jp) is 1 if p G db' and zero otherwise. Now since [n&'eE a dv\ = a Bb an< ^ smce 
ag b = (— l) ef, ( p ) we have that n(p) and e p (b) have to have the same parity proving that n(p) 
is even for p G Pq — db and odd for p G db. 
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Now lets take an element r G Ce- Since r(96') = 1 for all b' G E and E C E we have 
that: 

n = 1 

b'eE 

At the same time we can write: 

n = n = i 

b'<=E peP 
But since is even for p E Pq — db and odd for p G <% we have: 

n wp))" w = n t(p) = < B b) = i 

Thus we proved that if r G Ce then r(<%) = 1 and thus the new constraint is superfluous. 
This proves the lemma. 

Now using this lemma we see that the case where \Ce>\ = \\Ce\ and \T>e>\ = [De\ never 
happens and since in all the other cases our theorem holds then we proved our theorem. 



B.5 Summary 

We proved using theorem 2 that |X>e||Ce| > |^4| which means that \V\\C\ > \A\. Using this 
and the result of theorem 1: \C\\C\ = \A\ we see that \V\ > \C\ but we already know that 
since V C C we have that \V\ < \C\. The only possible solution to this is that \V\ = \C\ 
and since V C C we have that V = C. 

This is the result that we were looking for. Now we can say that: 

MCI aeV 

C[a\ = (15) 
1 a^V 

where V is the group formed by taking all possible products between the star configurations. 



C Appendix, Antiperiodic boundary conditions 

We have shown that vortices wrapped around the periodic boundary conditions have zero 
weight. However we point out how these configurations can instead be weight = 1 and the 
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above configurations weight = by a minor change in the formulation. Consider 

Z = J[dU]e(^ tr[Uidp)] ) - A^ e (^^^[^)]^E pe ,|tr [f /(a P )])_ (lg) 

where P is a co-plane of sign flipped terms in the action. This is known as antiperiodic or 
twisted boundary conditions in the literature. 

Antiperiodic boundary conditions amounts to nothing more than a change in the action. 
The derivation of the partition function in terms of SU(2)/Z2 and Z2 variables can be 
generalized. In our notation , Eqn.(|T6D becomes 

[du] e (pE P mu(d P) ]n P) ) ^ 



[dU] £ e 



Define: 



a' = a a <^> a = a a' 



where cxoQo) = — 1 for peP and +1 elsewhere. 
We can simplify Z, Eqn. (|i~7|) 

Z = J[dU] £ e (^^|trTO)lk'(,)) = j m ^ e (^E^|tr[^ P )] k »)_ 

This looks exactly like the partition function for the periodic boundary conditions except 
that r\a' G <j C instead of C. 
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